Fluorescent visualization of a spreading surfactant 
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The spreading of surfactants on thin films is an industrially and medically important 
phenomenon, but the dynamics are highly nonlinear and visualization of the surfactant dy- 
namics has been a long-standing experimental challenge. We perform the first quantitative, 
spatiotemporally- resolved measurements of the spreading of an insoluble surfactant on a thin 
fluid layer. During the spreading process, we directly observe both the radial height profile of 
the spreading droplet and the spatial distribution of the fluorescently-tagged surfactant. We 
find that the leading edge of spreading circular layer of surfactant forms a Marangoni ridge 
in the underlying fluid, with a trough trailing the ridge as expected. However, several novel 
features are observed using the fluorescence technique, including a peak in the surfactant 
concentration which trails the leading edge, and a flat, monolayer-scale spreading film which 
differs from concentration profiles predicted by current models. Both the Marangoni ridge 
and surfactant leading edge can be described to spread as R oc t s . We find spreading expo- 
nents, Sh ~ 0.30 and 6r ~ 0.22 for the ridge peak and surfactant leading edge, respectively, 
which are in good agreement with theoretical predictions of S = 1/4. In addition, we observe 
that the surfactant leading edge initially leads the peak of the Marangoni ridge, with the 
peak later catching up to the leading edge. 



I. INTRODUCTION 



The spreading of surfactant over a thin layer of fluid is of fundamental importance to diverse 
biomedical applications such as pulmonary surfactant replacement therapy [HE]) lipid tear layers 
in our eyes [3], and drug delivery mechanisms [4], as well as many industrial coating and draining 
flows. The dynamics of the spreading process are driven by highly nonlinear equations modeling 
the curvature of the fluid's free surface and gradients in the surfactant concentration [5H3]- Progress 
in assessing the efficacy of these theoretical models, recently reviewed by Craster and Matar \10\ 
111) , has been hampered by the difficulty of directly visualizing the location of the surfactant in 
experiments |12H14j . 

An experimentally convenient and commonly modeled geometry starts from a radially symmet- 
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ric droplet of surfactant which is allowed to spread over the underlying fluid [T2HT5] . The 
gradient in surfactant concentration at the leading edge of the droplet causes a gradient in surface 
tension. As a result, Marangoni forces pull outward on both the surfactant and the underlying 
fluid, causing the droplet to spread. During the spreading, a fluid ridge known as a Marangoni 
ridge forms near the leading edge of the surfactant. In models coupling lubrication theory for a 
thin fluid to the surface tension gradient provided by the surfactant concentration gradient, the 
Marangoni ridge and the peak or edge of the surfactant concentration profile have been associated 
with shocks or traveling wave solutions [SI H6HT8] . However, few experimental measurements have 
been made to quantify the relationship between the surfactant concentration and the fluid surface 
curvature, due to the difficulty of directly measuring the spatial distribution of surfactant molecules 

nam. 

Previous experiments have provided quantitative measurements of either the position of the 
surfactant leading edge or of the fluid surface profile, but not quantitative profiles of both for the 
same system. Bull et al. [13] measured the location of the leading edge of a spreading fluorescent 
surfactant on a glycerin layer, while fluorescent microspheres on the surface were used to infer 
the arrival of the surface disturbance. Dussaud et al. |14j . measured changes in the surface slope 
due to the spread of oleic acid on a glycerin-water layer via Moire topography and used talc 
particles on the surface to infer the arrival of the leading edge of the surfactant. In both studies, 
it was observed that the spreading surfactant lagged behind the surface disturbance, which we 
are calling the Marangoni ridge. In the experiments presented here, we provide simultaneous 
measurements of both the surface height profile and the surfactant concentration field by combining 
laser profilometry with further development of the fluorescence techniques introduced by |13] . This 
combination of techniques allows us to directly measure the spatiotemporal evolution of a droplet 
of a fluorescently-tagged surfactant lipid spreading on a thin layer of glycerin without introducing 
additional contaminants to the system. We find that, at later times, while the forward part of the 
surface disturbance can be said to be in front of the leading edge of the spreading surfactant, the 
surfactant leading edge coincides with the peak of the Marangoni ridge. 

We compare our experimental results with the models of Jensen and Grotberg [HI [JS] which 
predict a R ~ t 1 / 4 spreading law for the peak of the Marangoni ridge, similar to what was previously 
observed in experiments by Gaver and Grotberg [12J. (Note that this behavior is slower than that 
of droplets spreading on thick layers, for which the scaling is i 3 / 8 to i 3//4 , depending on the degree 
of surface contamination |19H21j or on very thin (< l//m) layers, for which the scaling is predicted 
and observed to be i 1/2 [3 [22].) In our experiments, we measure spreading exponents which 
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approximately agree with the predicted 1/4 values. However, the Marangoni ridge initially forms 
behind the leading edge of the surfactant, near the maximum in surfactant concentration, and 
only later catches up. Thus, the measured exponents are slightly higher for the position of the 
Marangoni ridge and slightly lower for the position of the leading edge. In the work of Dussaud 
et al. |14| . similar behavior was observed, namely a higher spreading exponent for the Marangoni 
ridge than for the surfactant leading edge. However, some care needs to be taken in connecting 
these two results, since the ridge position in |14j was tracked via the point of minimum surface 
slope rather than the peak. For a dispersive wave, as is the case here, the peak and the inflection 
point will become further apart over time. Similarly, the fluorospheres used to visualize the surface 
disturbance in Bull et al. [13] will be moved by the front of the wave associated with the surface 
compression disturbance, even before the peak of the Marangoni ridge arrives. By directly imaging 
both the full surface profile and the spatiotemporal evolution of the surfactant concentration, we 
avoid such difficulties of interpretation. 



The experimental cell (see Figure [Tp) consists of a machined aluminum well of radius L = 

111 mm, filled with 38 mL of 99.5% anhydrous glycerin; this results in a fluid layer of depth 
d = 0.98 ± 0.03 mm. Millimeter-scale fluid layers are chosen to allow for comparison with theories 
based on the lubrication approximation (d C i), while avoiding film rupture which would occur 
for thinner films |X2^ 114] . By using viscous (but still Newtonian) glycerin as the working fluid, we 
are able to achieve sufficiently slow dynamics at low Reynolds number to permit fluorescence mea- 
surements of the surfactant concentration. The glycerin is maintained at a constant temperature 
of 25 ± 0.5°C by circulating water from a temperature-controlled bath through the base of the alu- 
minum well; the viscosity of glycerin at this temperature is [i = 0.83 ± 0.03 Pa-s [23J. The surface 
tension of pure glycerin at this temperature is a$ = 63.4 dyne/cm [23]. The surface diffusivity of 
such molecules is small, approximately 10~ 4 cm 2 /s [25], and is frequently neglected in modeling 
the dynamics. Key material parameters and dimensionless ratios are summarized in Table [I] 

The spreading surfactant droplet is a fluorescently-tagged lipid from Avanti Polar Lipids, 1- 
palmitoyl- 2-12- [(7- nitro-2-l,3-benzoxadiazol-4-yl) amino] lauroyl-sn-glycero-3-phosphocholine (ab- 
breviated NBD-PC). This molecule is a two-tailed lipid with a hydrophilic (also glycerin-philic) 
head and hydrophobic 16-carbon and 12-carbon tails; the NBD fluorophore is attached to the 
shorter of the two tails which permits fluorescent visualization. NBD-PC is insoluble in glycerin, 



II. EXPERIMENTS 




4 




FIG. 1: (a) Schematic of initial configuration of experimental apparatus, side and top view, (b) Image of 
red laser profile of surface height and green fluorescence from surfactant at t = 2 sec, with green channel 
digitally enhanced rs 40 x and white dashed circle marking initial location of confining ring (exposure time 
1 second). The black rectangle denotes the area excluded from fluorescence measurements due to artifacts 
from the high intensity of the laser signal. See http://nile.physics.ncsu.edu/pub/movies/tliinfilms/ 



for a movie version of this figure, (c) Time-evolution of laser profile (image) with extracted profile H(r) 
(white line), for exposure time 0.1 sec. 

lowers the surface tension by as much as 30 dyne/cm |13| . has an absorption peak at 460 nm 
(excited by a broad-spectrum black light), and has an emission peak at 534 nm. 

The hydrophilic phosphocholine headgroup has a cross-section of approximately 50 A 2 |26j in an 
anhydrous medium; from this value, we estimate that ~ 0.1 mg of NBD-PC are sufficient to cover 
the full glycerin surface with a monolayer of molecules. This quantity corresponds to a surfactant 
monolayer concentration of Tq ~ 0.3 /jg/cm 2 and represents an upper limit on our monolayer 
concentration Tq because our liquid layer is not completely anhydrous. For phosphocholines in the 
presence of water, the cross-section would be no larger than 60 A 2 [26]. This critical concentration 
r c is important because the presence of a larger amount of surfactant will not further reduce 
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fluid thickness 


d 


0.98 mm 


lateral dimension 


L 


0.8 cm (ring radius) to 5 cm (final radius) 


aspect ratio e 


= d/L 


0.12 (initial) to 0.02 (final) 


fluid density 


P 


1.26 g/cm 3 


dynamic viscosity 


M 


0.83 Pa-s [23] 


bare glycerin surface tension 


^0 


63.4 dyne/cm [H] 


surfactant-contaminated surface tension 




35.5 dynes/cm [13] 


spreading parameter S = 


(7 — a m 


27.9 dynes/cm 


critical monolayer concentration 


r c 


« 0.3 /xg/cm 2 [13] 


surface diffusivity 


D 


10~ 4 cm 2 /s [25] 


Bond number (gravity/Marangoni forces) Bo - 


= pgh 2 /S 0.44 (0.16 for d = 0.59 mm) 


capillary/Marangoni forces e 


2 °o/S 


0.036 (initial) to 0.00091 (final) 


initial velocity 


V 


6 cm/s (maximum measured) 


Reynolds number Re 


— pvL/fi 0.73 (maximum measured) 



TABLE I: Summary of key dimensional and dimensionless parameters. 



the surface tension (i.e., Act ~ for V > T c ). This critical monolayer concentration is in good 
agreement with measurements by Bull et al. [13], where the equation of state ct(T) becomes flat 
for r > 0.3 //g/cm 2 . We conduct our experiments with concentrations which provide partial 
coverage of the full fluid surface. Most runs use m = 18 fiL of NBD-PC/chloroform solution which 
contains 18 /ig of NBD-PC. This amount corresponds to 1.3 x 10 16 molecules, which would provide 
a monolayer area with radius ~ 45 mm for a headgroup size of 50 A 2 . The experiments generally 
spread to approximately this area before we lose resolution in our fluorescence measurements. 

To perform surfactant-spreading experiments, we deposit the surfactant within a 16 mm diam- 
eter steel ring (4 mm tall) placed in the center of the glycerin layer. Prior to each experimental 
run the ring and the aluminum well were cleaned with Contrad soap and water then dried us- 
ing nitrogen gas. The surfactant is initially dissolved in chloroform; we allow at least 6 minutes 
of evaporation for each //L of the solution before lifting the ring to allow the lipid molecules to 
spread. The surfactant remains confined inside the ring during evaporation, and 18 //g are enough 
to provide approximately 30 monolayers within the ring. At the start of each experimental run, we 
lift the ring using a monofilament line connected by an overhead pulley to a slow motor to provide 
repeatable dynamics and minimize the formation of bubbles during release. 

Our laser profilometer is constructed from a 632.8 urn HeNe laser fitted with a line generator. 
We project a straight line of light across the glycerin surface at an angle of 12 ± 0.5° above the 
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horizontal. To create a thin line, the laser light passes through a 0.4 mm slit; this provides a 2 mm 
wide line at the surface of the glycerin. This line is refracted by the glycerin and reflected off the 
bottom surface of the aluminum well before being imaged by the camera overhead, providing a 
projection of the height profile of the glycerin surface (see Figure [lj>,c). We calibrate the laser 
deflection as a function of known displacements by noting the change in location of the laser line 
before and after adding a known thickness of glycerin to the well. Displacement of the laser line is 
a linear function of the surface displacement. 

We make simultaneous measurements of the fluid surface profile and the spatial surfactant 
distribution via laser profilometry and fluorescence, respectively. Using a color digital camera with 
a 525 nm high-pass filter, we can observe both the 632.8 nm laser profile and the 534 nm NBD-PC 
fluorescence while filtering out the UV light used for excitation. The height profile is recorded in 
the red channel of the camera and the surfactant fluorescence in the green channel, without the 
need to align two images. We utilize two different frame rates (exposure times) in collecting data 
in order to optimize either the fluorescence measurements or the fast dynamics. Using a 1 second 
exposure time, we are able to collect more emitted photons from the fluorescent surfactant and 
thus obtain improved spatial resolution. Using a 0.1 second exposure time, we are able to obtain 
greater resolution and earlier times in tracking the position of the Marangoni ridge. 

To determine the relationship between fluorescence intensities measured by the camera and 
surfactant concentration, we take pictures in which known quantities of NBD-PC are placed in 
20 mm diameter rings under normal experimental conditions. The choice of ring size is guided 
by a desire to minimize errors due to heterogeneity in NBD-PC concentration, yet still provide a 
large enough area for measurement. The intensities are averaged over a 12 mm diameter region, 
to avoid shadows at the edges of the ring. Over the range 0.14 to 8.8 /ig/cm 2 , we find that there 
is a linear relationship between intensity and NBD-PC concentration. We use this empirical linear 
relationship to convert fluorescence intensities to surfactant concentrations. Statistical errors in 
r are about ±0.05 /ig/cm 2 in the vicinity of the leading edge, and larger near the center of the 
droplet. 

III. FLUID AND SURFACTANT DYNAMICS 

In long-exposure images (1 second, as shown in Figure^), we can observe both the height profile 
and the surfactant concentration. Brighter green corresponds to greater surfactant concentration; 
the red line is a refracted laser line showing the height profile of the glycerin surface. Because 
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the fluorescence signal is weak, we use image-division to remove background heterogeneities and 
a median filter to eliminate the noisiest pixels. The Marangoni ridge, created by the spreading 
surfactant, is visible as the peak of the red laser line. The dashed circle indicates the location of 
the confining ring for the initial surfactant deposit. At early stages, as shown in Figure [TJ^-f, the 
spreading surfactant typically has a maximum concentration which trails the leading edge. This 
maximum is not present in the models, which predict a monotonic decrease [51 116 1 [271l29| . 

Figure [TJ3 shows examples of the raw laser line and the extracted height profile H(r) for short- 
exposure images (0.1 second). To determine the profile, we find the centroid of pixel intensity 
for each image column in the red channel of the image, apply a median filter with respect to 
neighboring columns, and smooth the resulting function using a local linear regression. Within 
each H{r) we fit a parabola to the Marangoni ridge to determine the location of the peak Rh ■ As 
the ridge moves outward from the center (r = 0) over time, the height of the peak decreases due 
to gravity and capillarity. A trough behind the ridge also develops but does not rupture [8j [H] 
the fluid layer for d > 0.5 mm. Eventually the peak moves out of the viewing area or becomes 
indistinguishable from the rest of the glycerin layer. As would be expected, the maximum height 
is observed to be less than twice the thickness of the fluid layer, the value predicted by models 
neglecting gravity and capillarity [51 1161 [27]. an d is consistent with models that include gravity [6]. 

To visualize the surfactant driving this motion, we use the long-exposure images (1 second) and 
radially average the intensity of the green channel for use in the empirically-measured conversion to 
concentration. Figure [2^,-f shows corresponding fluorescence concentration profiles T(r) and height 
profiles H(r) at representative times. Figure [2^, shows the first image after the ring lifts off the 
surface, with the surfactant still localized to the original ring location (and the shadow of the ring 
still visible). We set t = at the middle of the interval captured by the preceding image. As time 
progresses, a Marangoni ridge develops behind the leading edge of the surfactant and eventually 
catches up. Note that the fluorescence concentration within the spreading front has a plateau 
with a r ~ r c . Because the film has not ruptured, the surfactant does not become trapped at 
the point of deposition as in Dussaud et al. [14J, and the fluorescence concentration at the center 
continues to diminish as the surfactant spreads. Importantly, the measured shape of T(r) differs 
significantly from what has previously been predicted by various models of similar configurations 
[IH El HH [161 E2H25]. In models, T(r) falls off monotonically without the presence of either the 
plateau or the leading edge peak observed here. 

The distinct peak in surfactant concentration just behind the leading edge of T(r) is generated 
during the interval between the first two images (see Fig. ^b-b) and might therefore be associated 
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FIG. 2: (a-f) Comparison of fluorescence concentration (r(r), thin green, right axis) and height (H(r), 
thick red, left axis) profiles. The locations of Marangoni ridge peak Rh and surfactant leading edge Rr are 



marked with o. See http://nile.physics.ncsu.edu/pub/movies/thinfilms/ for a movie version of this 



figure, (g) Space-time plot of fluorescence concentration, overlaid with Marangoni ridge peak (solid red line), 
surfactant leading edge (solid green line), and initial ring radius (horizontal dashed black line) locations. 

with the detachment of the fluid from the ring. However, in control experiments in which the 
confining ring is lifted from a uniform layer of surfactant no such peak is formed. Because the 
surface slope is shallow (< 2°) in the vicinity of the peak, the increased intensity is unlikely to be 
due to an increased surface area arising from the tilt. The magnitude of the peak is not consistent 
between runs, and it also appears to be unstable since its magnitude decays in time (see Fig.[2]>f). 
Further work is necessary to understand how it arises and decays, and whether it can be reproduced 
in models. 

An additional unexpected feature is visible in Fig. [TJd: spatial heterogeneities in the surfactant 
concentration which appear as areas of brighter green between the initial ring location and the 
surfactant leading edge. These heterogeneities form during the initial deposit of the surfactant- 
chloroform mixture within the ring and do not dissipate due to the low surface diffusivity. The 
origin of the aggregation likely lies in the electrostatic attractions between the lipid molecules [30] . 
No characteristic length scale has been observed, but fingering instabilities in spreading surfactant 
fronts have been observed for flows down inclined planes [31] and for microscopic films [32H35J . 
Importantly, these heterogeneities likely have only a small effect on the spreading behavior since 
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FIG. 3: (a) Marangoni ridge location for eleven runs with m = 18 fig (thin solid lines) and Rn(t) oc t^ H ^ 
for the mean spreading exponent Sh measured for each run (dashed line), (b) Comparison of (Rait)} for 
several combinations of m and d: m — 9 fig and d — 0.98 mm (3 runs), 18 fig and d = 0.98 mm (11 runs 
shown in (a)), 36 fi and d — 0.98 mm (2 runs), and m = 18 fig and d — 0.59 mm (2 runs). Dashed black line 
shows comparison to spreading exponent S = 1/4. 

the spreading front is observed to remain circular at all times. This can be understood as arising 
from the fact that for T > T c , surface tension gradients are small. Additional work is needed to 
characterize this aggregation phenomenon, for which the ability to directly observe the surfactant 
location is a key technique. 

To visualize the spatial relationship between the Marangoni ridge and the surfactant concen- 
tration, Figure [2^ shows the positions of the peak of H(r) and the leading edge of T(r) with 
respect to the underlying fluorescence concentration field T(r, t). At early times (t < 6 seconds), 
the Marangoni ridge peak lags behind the surfactant leading edge. After about 6 seconds, the 
peak of the Marangoni ridge catches up to the leading edge of the surfactant, which corresponds 
to the behavior observed in [8]. The peak and leading edge move together for the remaining time 
we are able to track the position of the Marangoni ridge peak, which becomes too flat to detect 
after about 20 seconds. Therefore, no long-time comparison of the two positions is possible. 

To quantify the spreading of both the Marangoni ridge and the surfactant, we examine the 
scaling of the Marangoni ridge position Rait) and the leading edge position Rr(t). Figure [3^ 
shows Rn(t) for 11 experimental runs with m = 18 fig of NBD-PC, with each line representing a 
different experimental run. There is considerable scatter in the initial Rh(0), but the exponent 
is consistent between runs. We fit each curve to the form Rh oc t Sn and find a mean spreading 
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FIG. 4: (a) Position of the surfactant leading edge for three runs with m = 18 /ig. The dashed black line 
shows a comparison to spreading exponent <5r = 1/4. 

exponent of (5u) = 0.295 ± 0.03, represented by the slope of the dashed line. 

We examine the role of surfactant concentration by comparing runs with initial mass m of 9, 18 
and 36 fig. Figure [SJa shows (Rnit)) for each case, where the average is taken over the ensemble 
of runs. Interestingly, the spreading behaviors for m = 9 and 18 /ig are very similar, while the 
36 /ig has a larger prefactor, but still the same exponent 5h- It is possible that the faster rate 
is related to reservoir effects stemming from the additional supply of molecules, but in all cases 
the initial supply concentration exceeds one monolayer. For a thinner layer {d = 0.59 mm) which 
better-satisfies the conditions of the lubrication approximation without being so thin as to cause 
surface rupture, we observe a spreading exponent closer to 1/4. 

Figure [i] shows a similar plot for Rr(t) oc t Sr , using images with 1 second exposure times. We 
automatically located the leading edge by finding this peak in the intensity gradient of a smoothed 
profile, fitting a parabola to the peak, and taking the edge to be the position at which the parabola 
crosses zero. Our ability to detect the edge is limited at long times by the decreasing intensity of the 
signal. However, we do not see a change in behavior as the system approaches the estimated 45 mm 
radius at which the NBD-PC molecules become a monolayer. Figure [2]F supports the identification 
of this length scale, as the surfactant concentration profile T(r) has become nearly flat behind the 
leading edge at i?r = 40 mm. Three separate runs with m = 18 /ig yield (5r) ~ 0.22 ± 0.03, which 
is slower than observed for the Marangoni ridge position, and in approximate agreement with the 
1/4 scaling predicted by Jensen and Grotberg [8l [16] . 
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IV. DISCUSSION AND CONCLUSIONS 

The combined techniques of laser profilometry and fluorescent surfactants have allowed us to 
simultaneously monitor both the surface deformation of a glycerin layer and the spreading of 
surfactant molecules across the surface. Releasing a mass of surfactant from within a confining ring 
allows us to compare our experimental observations with models and simulations for the spreading 
of an axisymmetric drop as in [HI US] which predict R oc i 1 / 4 . We find that the measured exponents 
(^r) = 0.22 ± 0.03 and (5h) = 0.295 ± 0.03 are both in reasonable agreement with the predicted 
exponent. It is important to note that the exponent 1/4 is in fact only expected under a set of 
highly restrictive conditions which neglect gravity, neglect capillarity, and utilize the lubrication 
approximation (e <C 1). Furthermore, the lifting of the ring creates a disturbance with its own 
length scale: it is therefore likely that no single exponent fully characterizes the spreading. As 
observed in the models, the peak of the Marangoni ridge and the location of the surfactant leading 
edge coincide (although only after an initial transient during which the ridge lags). In addition, 
the fluid profile is observed, as in [6] , to have a capillary trough which trails the Marangoni ridge. 
The fluorescence imaging indicates that this trough is associated with a reservoir for the spreading 
surfactant. 

Interestingly, the spreading exponent 5h (measured for the peak of the Marangoni ridge) is 
higher than 5 = 1/4 but <5r (obtained from fluorescence measurements of the position of the leading 
edge of the surfactant) is lower. While the range of times over which we can take measurements 
of Rh and Rp are overlapping, early-time data is only available for Rh and late-time data is only 
available for Rp. Since Rr is measured over both longer and later times than Rh, <5r would be 
expected to be closer to the asymptotic result. The higher-than-expected value of 5h may simply 
reflect the lingering effects of the initial disturbance provided by the ring lifting and then detaching 
from the fluid surface. Thus, its possible that an (unobserved) crossover to 5r = Sh occurs for 
later times than were measured. 

In addition, several other discrepancies from the idealized model [8] [16] might act separately on 
the two positions Rh and Rt- Capillarity and gravity do not act on the fluid and surfactant in the 
same way. Second, while the ratio of film thickness to lateral extent is small, three-dimensional 
effects (such as circulation under the free surface) are likely still present and may affect the spread- 
ing rate of the fluid (and to a lesser extent, the surfactant). Indeed, Gaver and Grotberg |12j 
calculated that there is a re-circulation of particle-trajectories within the bulk of a fluid, but it 
remains unclear what consequence this has on the Marangoni ridge motion. Further efforts to 
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probe the validity of the lubrication approximation in this regime are necessary. 

To understand the low value of <5r, we should consider the surfactant itself. In our initial 
configuration, the confining ring contains multiple layers of surfactant. These multiple layers may 
act as a reservoir of surfactant and thus might provide an effective flux of surfactant. However, the 
predicted effect [8] is 5 = (1 + ct)/4 for a constant flux of surfactant supplied at a rate of t a , which 
would necessarily lead to a larger value for 5r, rather than the smaller value observed. In addition, 
we do not observe a change in 5h ~ 0.3 even when we change the initial number of monolayers 
via initial surfactant mass m. A second effect is the presence of a non-zero surface diffusivity [25j. 
which is predicted to increase 5 by amounts similar to the discrepancies observed here [B]. However, 
this would shift Sr in the opposite direction from what is observed. 

Significantly, two key discrepancies with existing models [6j 00 [131 03 [161 [27] ar i se m examining 
the shape of the surfactant concentration T(r): neither the presence of a plateau behind the leading 
edge nor the peak immediately behind the leading edge are predicted. Instead, the model T(r) 
falls off monotonically. Interestingly, the plateau has a near- monolayer concentration (r ~ T c ). At 
early times, the region behind the leading edge has T > T c , and therefore only small surface-tension 
gradients would act to smooth the peak in T. Thus, there is no observable disturbance to H(r) 
due to peak in T(r). Further modeling work will be required in order to understand how both the 
extended plateau and the peak arise and ultimately decay. 
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anonymous reviewers. Keith Weninger provided a spectrophotometer for use in measuring the 
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